
capture log close
log using $data/EDS_08_robustness_intensity, text replace

clear
clear matrix
clear mata
set mem 2g
set matsize 800
set more off
set maxvar 20000
set trace off
set varabbrev off

#d;
	  global cov "_Isu_schule_1 _Isu_schule_2 _Isu_schule_3 
	_Isu_ausbila1 _Isu_ausbila2
	_Isu_kat_al_2 _Isu_kat_al_3 _Isu_kat_al_4 
	_Ichild_1 _Ichild_2
	_Isu_eintri_1 _Isu_eintri_2 _Isu_eintri_3 _Isu_eintri_4 _Isu_eintri_5
	_Isu_eintri_6 _Isu_eintri_7 _Isu_eintri_8 _Isu_eintri_9 _Isu_eintri_10 _Isu_eintri_11
	_Ibefstat_1 _Ibefstat_2 _Ibefstat_3 _Ibefstat_4
	_Idauloint__8 _Idauloint__9 _Idauloint__10 _Idauloint__11 _Idauloint__12 
	_Idauloint__13 _Idauloint__14
	female east_germany
    german migration_background 
    married  
	benefit_status benefit_level log_benefit_level
    lifetime_ue_age lifetime_empl_age last_wage last_wage_mis 
    region_ba_aloq_zivall gdp_real_pc region_ba_vacancies 
	sector_manufacturing sector_service  log_mobirate_avg5
	openness_st conscient_st extraversion_st neuroticism_st loc_st
	write_ger speak_ger write_eng speak_eng
	homeowner   
	   ";		

	   

		global outc"
		   search_dist50
		   log_distance_km 
	       applications_total
		   applications_distant
		   applications_local
		   
		   
		   ";
	  #d cr
	  
	  
	  use $data/data_pair, clear
	
	  
	  	
		gen intensity = exp(log_intensity_alt)
		
		gen intensity_high=0
		forvalues x=1/487 {
			qui sum log_intensity_alt if county_pair_`x'==1
			replace intensity_high=1 if log_intensity_alt>r(mean) & county_pair_`x'==1
		}
		
		gen intensity_high25=0
		forvalues x=1/487 {
		qui sum log_intensity_alt if county_pair_`x'==1 
		qui replace intensity_high25=1 if log_intensity_alt>(r(mean)+0.25) & county_pair_`x'==1 
		}
		
		
		
		
		
		gen log_intensity_inc = log(intensity / hheink) 
		
	    foreach type in  $outc {
	 	
			
		
		reg `type'  log_intensity_alt  $cov  county_pair_*  [pw=1/fweight_reg], cluster(region)
	    matrix b=e(b)'
		matrix n=e(N)
		matrix v==e(V)
		matrix se=vecdiag(cholesky(diag(vecdiag(v[1..1,1..1]))))'
		matrix d0_`type'=b[1..1,1..1],se
		
		
		qui reg `type'  log_intensity_inc $cov  county_pair_*  [pw=1/fweight_reg], cluster(region)
	    matrix b=e(b)'
		matrix n=e(N)
		matrix v==e(V)
		matrix se=vecdiag(cholesky(diag(vecdiag(v[1..1,1..1]))))'
		matrix d1_`type'=b[1..1,1..1],se
				
			
		qui reg `type'  log_intensity_exits $cov  county_pair_*  [pw=1/fweight_reg], cluster(region)
	    matrix b=e(b)'
		matrix n=e(N)
		matrix v==e(V)
		matrix se=vecdiag(cholesky(diag(vecdiag(v[1..1,1..1]))))'
		matrix d2_`type'=b[1..1,1..1],se
			
			
		qui reg `type'  intensity $cov  county_pair_*  [pw=1/fweight_reg], cluster(region)
	    matrix b=e(b)'
		matrix n=e(N)
		matrix v==e(V)
		matrix se=vecdiag(cholesky(diag(vecdiag(v[1..1,1..1]))))'
		matrix d3_`type'=b[1..1,1..1],se
		
		qui reg `type'  intensity_high $cov  county_pair_*  [pw=1/fweight_reg], cluster(region)
	    matrix b=e(b)'
		matrix n=e(N)
		matrix v==e(V)
		matrix se=vecdiag(cholesky(diag(vecdiag(v[1..1,1..1]))))'
		matrix d4_`type'=b[1..1,1..1],se
				
		qui reg `type'  intensity_high25 $cov  county_pair_*  [pw=1/fweight_reg], cluster(region)
	    matrix b=e(b)'
		matrix n=e(N)
		matrix v==e(V)
		matrix se=vecdiag(cholesky(diag(vecdiag(v[1..1,1..1]))))'
		matrix d5_`type'=b[1..1,1..1],se
				
		
		matrix d1_`type'=d0_`type',d1_`type',d2_`type'
		
		matrix d2_`type'=d3_`type',d4_`type',d5_`type'
		
		
		
		
	
}/*type*/
		#delimit;
		matrix d1=
		d1_search_dist50 \
		d1_log_distance_km \ 
	    d1_applications_total \
		d1_applications_distant \
		d1_applications_local
		;
	
		matrix rownames d1 = 
		"Applied to distant" 
		"Search radius" 
		"#applications total" 
		"#applications distant" 
		"#applications local" 
		;
		#delimit cr
		
		#delimit;
		matrix d2=
		d2_search_dist50 \
		d2_log_distance_km \ 
	    d2_applications_total \
		d2_applications_distant \
		d2_applications_local
		;
	
		matrix rownames d2 = 
		"Applied to distant" 
		"Search radius"  
		"#applications total" 
		"#applications distant" 
		"#applications local" 
		;
		#delimit cr

		matrix colname d1 = "Baseline" "SE" "by household income" "SE" "by exits from UE"  "SE"
		estout matrix(d1, fmt(4)) using $data/Table_A6_survey_upper.tex, style(tex) replace
		matrix colname d2 = "Level effect" "SE" "Indicator" "SE" "Indicator 25%+" "SE"
		estout matrix(d2, fmt(4)) using $data/Table_A6_survey_lower.tex, style(tex) replace
	
	
log close	

